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Abstract 



These lectures give an overview of the uses of effective field theories in de- 
scribing gravitational radiation sources for LIGO or LISA. The first lecture reviews 
some of the standard ideas of effective field theory (decoupling, matching, power 
counting) mostly in the context of a simple toy model. The second lecture sets up 
the problem of calculating gravitational wave emission from non-relativistic binary 
stars by constructing a tower of effective theories that separately describe each scale 
in the problem: the internal size of each binary constituent, the orbital separation, 
and the wavelength of radiated gravitons. 



1. Lecture I 



1.1. Introduction and motivation 

These lectures describe the uses of effective field theory (EFT) methods to solve 
problems in gravitational wave physics. Many of the signals relevant to gravita- 
tional wave experiments such as LIGO (TJ, VIRGO ||2) and the planned LISA |[3l 
correspond to astrophysical sources whose evolution involve a number of dis- 
tinctly separated length scales. In order to compute signal templates that capture 
the physics accurately, it is important to systematically account for effects arising 
at all these different scales. It is exactly this sort of problem that is best treated 
by EFT methods, analogous to the EFTs constructed to unravel multiple scale 
problems in high energy physics or condensed matter. 

For the sake of concreteness, I will focus in these lectures on the EFT formu- 
lation of the slow "inspiral" phase of compact binary stars (that is, with neutron 
star (NS) or black hole (BH) constituents). A more comprehensive review of 
gravitational wave sources and phenomenology can be found in the lectures by 
A. Buonnano at this school (also see, e.g., ref. (4)). The inspiral phase plays an 
important role in gravitational wave physics, and corresponds to the period in the 
evolution of the binary in which the system is non-relativistic, the bound orbit 
slowly decaying due to the emission of gravitational radiation. 

An appealing feature of the binary inspiral phase is that, theoretically, it is a 
very tractable problem. In principle, the dynamics is calculable as a perturbative 
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expansion of the Einstein equations in the parameter v -C 1, a typical three- 
velocityQ To get a handle on the scales involved, it is worthwhile to calculate 
some of the features of the binary inspiral at the order of magnitude level. Let's 
consider inspiral events seen by LIGO, which for illustration we will assume 
operates in the frequency band between 10 Hz and 1 kHz. At such frequencies, 
the relevant inspiral sources correspond to neutron star binaries, with m ns ~ 
1.4m0, (m Q is the solar mass) or perhaps black holes with typical mass itibh ~ 
IOitiq. Since the time evolution is non-relativistic, it is to a good approximation 
treated in terms of Newtonian gravity. In particular, the orbital parameters are 
related by 

, 2 ~^^, (l.D 
r ir 

where r is the orbital radius, and we have introduced the Schwarzschild radius 
r s = 2Gmm. The frequency of the radiation emitted during this phase is of 
order the orbital frequency 2irv ~ v/r, so for a binary inspiral that scans the 
LIGO band, the typical orbital radius is measured in kilometers 

/ \ V3 / \ 1/3 

/ 777 \ / 777, \ 

r(10 Hz) ~ 300 km -> r(lkHz)-14km ,(1.2) 

V m ©/ V m ©/ 

which is therefore well separated from the size of the compact objects, r s ~ 
1 kmm/mQ. The expansion parameter v also evolves as the binary sweeps the 
detector frequency band: 

/ N 1/3 / v 1/3 

/ 777 \ / 777 \ 

w(10Hz)-0.06 -> u(lkHz)~0.3 , (1.3) 

V m 0/ \m&J 

indicating that the velocity is a good expansion parameter up to the last few cycles 
of orbit seen by the detector. As v — > 1, the inspiral phase ends, and the evolution 
must be treated by numerical simulations. 

During the inspiral phase, the bound orbit is unstable to the emission of gravi- 
tational radiation. The power emitted in gravitational waves is well approximated 
by the quadrupole radiation formula. Taking the orbit to be circular, the power 
radiated is 



dE _ 32 j 1(l 
dt ~ 5 N 



G^v w . (1.4) 



1 Even though we will be doing classical physics for most of these lectures, I use particle physics 
units c = h = 1 throughout. 
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The mechanical energy of the binary is just E = —^mv 2 to leading order in v, 
so by conservation of energy 

| (-!„,„*) . _f G -v», (i,, 

we obtain an estimate for the duration of the inspiral event seen by LIGO 



1 1 

vf in 



, -8/3 



5min. , (1.6) 

m 0/ 



and the number of orbital cycles as the signal sweeps the detector band 



N ~ £ U {t)dt = i 



1 1 



4 x 10 4 ( — ) radians, (1.7) 



with the orbital angular frequency. 

Eqs. ( 11.2b . ( 11.3b . (11.6b give an estimate for the typical length, velocity, and 
duration of an inspiral event that sweeps the LIGO frequency band. Inspiral 
events seen by LISA follow the same dynamics. However, because LISA will 
operate in a frequency range complementary to LIGO, 10 -5 Hz < v < 1 Hz, 
the sources correspond to objects of much large mass (e.g., BH/BH binaries with 
rriBH ~ 1O 5 ~ 8 TO orBH/NS systems with m B H ~ 10 5 ~ 7 ). 

For either LIGO or LISA, the number of orbital cycles spent in the detector 
frequency band is quite large, see e.g. Eq. ( 11.71 ). Thus even a slight deviation 
between theoretical calculations of the gravitational wave phase (the waveform 
"templates") and the data will become amplified over the large number of cy- 
cles of evolution. Consequently, inspiral wave signals carry detailed information 
about gravitational dynamics. In fact, it has been determined that LIGO will be 
sensitive to corrections that are order v 6 in the velocity expansion beyond the 
leading order quadrupole radiation predictions 0. 

The procedure for computing corrections to the motion of the binary system 
in the non-relativistic limit by an iterated expansion of the Einstein equations is 
called the post-Newtonian expansion of general relativity [6|. What makes the 
calculations difficult (and interesting) is that there is a proliferation of physically 
relevant effects occurring at different length scales, 

r s = size of compact objects, 
r = orbital radius, 
A = wavelength of emitted radiation, 
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all controlled by the same expansion parameter v: 




(1.8) 



where the first estimate follows from Kepler's law and the second from the multi- 
pole expansion of the radiation field coupled to non-relativistic sources. Thus to 
do systematic calculations at high orders in v, it is necessary to deal with physics 
at many different scales. A convenient way of doing this is to construct a chain 
of EFTs that capture the relevant physics at each scale separately. 

The goal of these lectures is to recast the problem of computing gravitational 
wave observables for non-relativistic sources in terms of EFTs. In the next sec- 
tion we briefly review the EFT logic, using a standard particle physics example, 
the low energy dynamics of Goldstone bosons in a theory with spontaneously 
broken global symmetry, to illustrate the main ideas. In sec. [2] we turn to the EFT 
formulation of the binary inspiral problem. The presentation follows [7], which 
is based on similar EFTs for non-relativistic bound states in QED and QCD 00. 
However, the focus here is not on detailed calculations, but in describing how to 
"integrate out" physics at each of the scales r s , and r to obtain a theory with well 
defined rules for calculating gravitational wave observables (at the scale A) as an 
expansion in v. Possible extensions of the ideas presented here to other problems 
in gravitational wave physics are discussed in the conclusions, sec. [3] 

1.2. Effective field theories: a review 

EFTs are indispensable for treating problems that simultaneously involve two or 
more widely separated scales. As a typical application, suppose one is interested 
in calculating the effects of some sort of short distance physics, characterized by 
a scale A, on the dynamics at a low energy scale ui <C A. In the EFT description 
of such a problem, the effects of A on the low energy physics become simple, 
making it possible to construct a systematic expansion in the ratio uj/A <C 1. 
Here we briefly review the basic ideas that go into the construction and applica- 
tions of EFTs. More detailed reviews can be found in refs. [|9l ITOl ITTl IT2l . 

The key insight that makes EFTs possible is the following: consider for in- 
stance a field theory with light degrees of freedom collectively denoted by cf> (e.g, 
a set of massless fields) and some heavy fields $ with mass near the UV scale A. 
The interactions of these fields are described by some action functional S[<j>, <£>]. 
Suppose also that we are only interested in describing experimental observables 
involving only the light modes <fi. Then it makes sense to integrate out the modes 
$ to obtain an effective action that describes the interactions of the light fields 
among each other 




(1-9) 
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It turns out that in general, S e g [tf>] can be expressed as a local functional of the 
fields 4>{x), 



for (in general an infinite number) local operators Oi(x). The coefficients c, are 
usually referred to as "Wilson coefficients". If Oi(x) has mass dimension A», 
then the Wilson coefficients, evaluated at a renormalization point fi of order A, 
scale as powers of A 



with on ~ 0(1 ). From this observation we conclude that the short distance 
physics can have two types of effects on the dynamics at energies u> <C A: 

• Renormalization of the coefficients of operators with mass dimension A < 4. 

• Generation of an infinite tower of irrelevant (i.e. A > 4) operators with 
coefficients scaling as in Eq. 11.11) . 

This result regarding the structure of the low energy dynamics encoded in 
S e ff [<f>] is usually referred to as decoupling. This concept has its roots in the work 
of K. Wilson on the renormalization group lfl3ll . Decoupling as used by practic- 
ing field theorists was first made rigorous in [14]. The idea is useful because it 
states that the dependence of the low energy physics on the scale A is extremely 
simple. All UV dependence appears directly in the coefficients of the effective 
Lagrangian, and determining the A dependence of an observable follows from 
power counting (essentially a generalized form of dimensional analysis). Effec- 
tive Lagrangians are typically used in one of two ways: 

1. The "full theory" S[(j>, $] is known: In this case, integrating out the heavy 
modes gives a simple way of systematically analyzing the effects of the heavy 
physics on low energy observables. Because only the low energy scale ap- 
pears explicitly in the Feynman diagrams of the EFT, amplitudes are easier to 
calculate and to power count than in the full theory. 

Examples: Integrating out the W, Z bosons from the SU(2)l x U(1)y Elec- 
troweak Lagrangian at energies E <C mw,z results in the Fermi theory of 
weak decays plus corrections suppressed by powers of E 2 /rriyy z . It is easier 
to analyze electromagnetic or QCD corrections to weak decays in the four- 
Fermi theory than in the full Electroweak Lagrangian, as the graphs of Fig.Q] 
indicate. See, e.g., ref. lfT31 . Another example is the use of EFTs to calculate 
heavy particle threshold corrections to low energy gauge couplings lfl6l . This 
has applications, for instance, in Grand Unified Theories and in QCD. 




(1.10) 



Ci(fi = A) 



AAi-4' 



(l.H) 
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2. The full theory is unknown (or known but strongly coupled): Whatever the 
physics at the scale A is, by decoupling it must manifest itself at low ener- 
gies as an effective Lagrangian of the form Eq. (1 1 . 10b . If the symmetries (eg. 
Poincare, gauge, global) that survive at low energies are known, then the op- 
erators Oi(x) appearing in S e ff[4>] must respect those symmetries. Thus by 
writing down an effective Lagrangian containing the most general set of op- 
erators consistent with the symmetries, we are necessarily accounting for the 
UV physics in a completely model independent way. 

Examples: The QCD chiral Lagrangian below the scale A x sb of SU(3)l x 
SU(3) r -> SU(3) V chiral symmetry breaking El US). Here the full the- 
ory, QCD, is known, but because A x sb is of order the scale Aqcd ~ 1 GeV 
where the QCD coupling is strong, it is impossible to perform the functional 
integral in Eq. ( 12.161 ) analytically. Another example is general relativity be- 
low the scale mpi ~ 10 19 GeV. This theory can be used to calculate, e.g., 
graviton-graviton scattering at energies E <C mpi. Above those energies, 
however, scattering amplitudes calculated in general relativity start violating 
unitarity bounds, and the effective field theory necessarily breaks down. Thus 
general relativity is an effective Lagrangian for quantum gravity below the 
strong coupling scale mpi. The EFT interpretation of general relativity is 
reviewed in more detail in refs. Ifl9ll20l 

Finally, it is believed that the Standard Model itself is an effective field theory 
below scales of order A = 1 TeV or so (see lectures by H. Murayama at this 
school). This scale manifests itself indirectly, in the form of SU (2)l x U(1)y 
gauge invariant operators of dimension A > 4 constructed from Standard 
Model fields ETI . certain linear combinations of which have been constrained 
experimentally using collider data from the LEP experiments at CERN and 
from SLD at SLAC (see lf22l for a recent analysis of precision electroweak 
constraints using effective Lagrangians). If there is indeed new physics at the 
TeV scale, it will be seen directly, at the CERN LHC which is due to come on 
line in the next few years. 

In either of these two classes of examples, integrating out the heavy physics 
as in Eq. ( 12.161 ) results in a an effective Lagrangian that contains, in general, an 
infinite number of operators Oi{x). However, because 

1. An operator 0( x ) with [O] = Ae> > 4 contributes to an observable at relative 
order 



2. A given observable can only be determined up to a finite experimental resolu- 
tion e <C 1, 




(1.12) 
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one may typically truncate the series in Eq. ( II. 10b after a finite number of opera- 
tors, those with mass dimension A < N + 4, where 

e = expt. error = y— J . (1-13) 

Therefore S e f / [<fi) is predictive as long as there are more observables than oper- 
ators with mass dimension A < N + 4. 

To make these general remarks more concrete, it is worth studying a toy exam- 
ple in some detail. We'll consider the quantum field theory of a single complex 
scalar field 4>(x) with Lagrangian (the "full theory") 

C = \d^\ 2 - V ((/>), (1.14) 

where 

n<«^(H 2 -«72) 2 . (i.i5) 

This theory is invariant under the global U(l) symmetry 

U(l) : 4>{x) -> e la (j){x), (1.16) 
for some constant phase a, as well as a discrete charge conjugation symmetry 

C : <t>{x) -> 4>*{x). (1.17) 

Classically, the ground state of this theory is determined by the minimum of 
V(4>). Because of the 17(1) symmetry the ground state is degenerate, and the 
vacuum manifold is the circle \(j>\ 2 — v 2 /2 in the complex <j> plane. To study 
fluctuations about the vacuum, we expand the fields about any of these (equiva- 
lent) vacuua, for instance the point 

W = -T=- ( L18 > 

Expanding about this point spontaneously breaks the U(l) symmetry, resulting 
in one Goldstone boson. It is convenient to write the original field <p(x) as 

<j){x) = ^={v + p{x))e™ {x)/v . (1.19) 
v 2 

Under the original symmetries, the new fields transform as 

U(l) ■ { p(x) p(x) > (1 20) 
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and 

r . J p( x ) -> P( x )i n ?n 

C ■ { TT(X) -> -7T(l). (L21) 

The Lagrangian in terms of the new fields is given by 

£ = \{d»P? + ~ (l + ~)\d^f ~^(vp + p 2 /2) 2 , (1-22) 

so the spectrum of excitations about the vacuum consists of a "modulus" field 
p(x) with tree level mass given by m 2 p — \v 2 and the Goldstone boson tt(x) 
with = 0. (We take A <C 1 so that a perturbative treatment is valid). 

Suppose that we are interested in working out the predictions of this theory at 
energies u> -C m p . At such scales, only the dynamics of the massless field tt(x) 
is non-trivial, and following our general discussion it is convenient to write an 
effective Lagrangian. The general structure of this effective Lagrangian is dic- 
tated by the symmetries of the original theory. In particular the U(l), realized 
non-linearly as the shift symmetry tt(x)/v — > tt(x)/v + a, restricts the effective 
Lagrangian to be a function of d^ir only. Therefore the symmetries of the low en- 
ergy theory alone explain a number of consequences of the low energy dynamics, 
for instance 

• The field ir(x) must be massless, since a mass term would break the shift 
symmetry. This is just the statement of Goldstone's theorem in the context of the 
low energy EFT. 

• The field ir(x) is derivatively self-coupled. This implies in particular that 
scattering amplitudes are "soft", vanishing as powers of the typical energy u> in 
the limit uj — > 0. 

In addition to the constraints from the non-linearly realized U(l), there is also a 
constraint from the charge conjugation symmetry tt — > — ir which says that the 
EFT must be even in n. Thus the effective Lagrangian must be of the general 
form 

Ceft = \{d,nf + ^(^tt^tt) 2 + •• • , (1.23) 

where only the leading operators with two and four 7r's have been displayed. In 
this equation A ~ m p and c$ is some dimensionless constant of order A. 

7.2.7. Matching 

Since we know the full theory, Eq. (11.14b . it is possible to explicitly perform the 
functional integral of Eq. (11. 9\ in order to obtain the EFT parameters like A and 
c§ in terms of the couplings A, v appearing in the full theory Lagrangian. Rather 
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than calculating this integral, the low energy parameters are fixed in practice 
through a procedure referred to as matching. 

To perform a matching calculation, one simply calculates some observable, 
for instance a scattering amplitude among the light particles, in two ways. First 
one calculates the amplitude in the full theory, expanding the result in powers of 
u>/A. One then calculates the same quantity in the effective field theory, adjusting 
the EFT parameters in order to reproduce the full theory result. 

As an example, consider irir — > irir scattering in our toy model. In the full 
theory, the amplitude to leading order in A is given by the diagrams of Fig. [2] 
Reading off the Feynman rules from the full theory Lagrangian, Eq. ( 11.141 ). one 
finds 



i-Afuii = 





i 


' 2 H t" 


— Ki • k 2 




— fc 3 • fc 4 




V 


s — m 2 p 


V 



crossings, 



(1.24) 



where k\, ki are the initial (incoming) momenta, and k$, k^ are the final (out- 
going) momenta. Introducing the usual Mandelstam variables s = (hi + ki) 2 , 
t = (ki — ka.) , u = (ki — kA 



iAfuii 



fc 3j r, u 
u 

v 2 m 2 p 



(fci- 



2 , and expanding in s,t,u m 2 , this becomes 
0(s 3 /m 6 p ) (1.25) 



UN 2 

2 



In the effective field theory, the amplitude arises from the leading dimension eight 
interaction in Eq. dl.23l ). 



lAEFT = TV 



A 4 



(I) 1 



(1.26) 



Thus taking A = m p , the matching condition for the coefficient c 8 is given by 



c s (fi = A) = 4A + C(A 2 ). 



(1.27) 



7.2.2. Corrections to the matching coefficients 

In general, the coefficients in the effective Lagrangian are calculable as a series 
expansion in the parameters of the full theory. In our example, the coefficients 
receive corrections at all orders in the parameter A <C 1, of which the tree level 
matching calculation in the above example gives only the first term. It is instruc- 
tive to see how one goes about computing corrections to the Wilson coefficients, 
as this illustrates certain general properties of effective field theories. 

In the example above, consider for instance the coefficient of the kinetic term 
for the 7r field in the low energy EFT, 

Ceft = \z{d^f + ■ ■ ■ . (1.28) 
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g 




(a) (b) 



Fig. 1. QCD correction to a typical four-Fermi process in the Standard Model, calculated in (a) 
the full theory with propagating W bosons, and (b) the effective Fermi theory of weak interactions. 
Graph (b) reproduces graph (a) up to corrections suppressed by powers of E 2 /m 2 v <C 1. 



1\ / 3 

\ / 
\ / 
\ / 
\ / 
\ / 

\ ^ + crossings 

/ \ 
/ \ 
/ \ 
/ \ 
/ \ 
/ \ 

1/ \4 



Fig. 2. Leading order contribution to nir — > nir scattering in the full theory. The intermediate double 
line corresponds to the p field propagator. 
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(a) 




(c) 

Fig. 3. Feynman graphs contributing to the n field self-energy in the full theory with a propagating p 
field. 



To leading order Z = 1. Corrections of order A are obtained by comparing the 
one-loop 7r field propagator in the full and the effective theories, adjusting Z 
so that the two calculations agree. Rather than calculate the full propagator, it 
is enough to consider the "one light-particle irreducible" two-point function. In 
the full theory this is the sum of graphs that cannot be made disconnected by 
cutting one tt field propagator. In the EFT, it is just the usual 1PI graphs. The 
relevant diagrams are shown in the full and effective theories in Fig. [3] and Fig. [4] 
respectively. In the full theory we find at one-loop 



Fig.EIa) 



iX 
16tt 2 




+ 0(k 4 /m 2 p ), 



(1.29) 



Fig. Kb) = - 



iX 
AM 2 



7 + 1 - In 

e 



(1.30) 



where the graphs have been regulated by dimensional regularization in d = 4— 2e 
dimensions. The third diagram is 



Fig.Sc) 



2iXk 2 



d d q 
(2n) d 



0, 



(1-31) 



where we have made use of the fact that in dimensional regularization 



lim 



d d q 1 



= lim 



T(n-d/2) 



^oj (2n) d {q 2 ) n rwo (47r) d / 2 T(n) 



(0) 



d/2-n 



(1.32) 
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i \ 
i i 
\ / 

_ v _w-_ 

(a) (b) 

Fig. 4. Feynman graphs contributing to the n field self-energy in the low energy EFT. Graph (b), 
including one insertion of the coupling eg vanishes in dimensional regularization. 



The last diagram is given by 



Fig.[3d)=-|^fc 2 



1 

- - 7" 
e 



1 - In 



(1.33) 



Note that the 1/e divergences in Fig.[3ja) and Fig.[3tb) cancel each other. This 
is equivalent to the statement that the full theory does not have wavefunction 
renormalization at one-loop order (it is just </> 4 after all). On the other hand, 
the tadpole graph in Fig.[3jd), which corresponds to a radiative correction to the 
VEV of the original field <fi requires a tree graph (not shown in the figure) with 
a counterterm insertion to regulate the logarithmic divergence. After subtracting 
this 1/e pole using the counterterm contribution, one finds the following result 
for the renormalized tt field self-energy in the full theory at one-loop (working in 
the M S scheme): 



- ' - lX k 2 



16tt 2 




G(fc 4 /m 2 ). 



(1.34) 



In the EFT, the self-energy of the Goldstone boson is given by the graphs of 
Fig- El These are: 



Fig.Ha) = i(Z-l)fc 2 , 
and 

Fig. lib) = 0, 

by Eq. ( 11.321 ). Setting the full and EFT results equal, we find 



Z = 1 



A 



16tt 2 



?-31n 



0(X 2 



(1.35) 



(1.36) 



(1-37) 



The one-loop matching calculation presented here serves to illustrate two gen- 
eral features of matching calculations: 



Les Houch.es Lectures on Effective Field Theories and Gravitational Radiation 1 5 



• Integrating out heavy modes not only generates new irrelevant (dimension 
A > 4) operators in the effective theory. It also renormalizes the coefficients of 
the operators with A < 4. 

• Loop graphs in matching calculations typically contain logarithms of m p //x. 
In order to avoid possible large logarithms that could render perturbation theory 
invalid, one must choose a matching scale fi that is of the same order as the 
masses of the fields that are being integrated out. 

This last point in particular implies that it is the Wilson coefficients Cj (fi ~ A) 
that exhibit simple scaling in powers of the UV scale A. Given the RG equations 
for the full theory coupling constants, one may also use any other /i > A as the 
matching scale. Possible large logarithms that may arise from this choice are 
then resummed by RG running. For example, using 



^ = 3A(m) + 



and the one-loop RG equation for the full theory coupling constant A, 

one finds the relation 

Z( M ) = Z(A) + jjln^j. (1.40) 

Thus in this theory, using the "wrong" matching scale has a relatively mild effect, 
although this is not generally true in other theories, e.g. QCD. In any case, it is 
simplest in practice to just match directly at /i ~ A. 

Likewise, to calculate observables at low energy scale uj -C A, it is better to 
evaluate loop graphs in the EFT at a renormalization point fi ~ us, indicating that 
in the EFT one should use the couplings c;(/i ~ uj). These can be obtained in 
terms of the coefficients Cj (fi ~ A) obtained through matching by renormaliza- 
tion group (RG) evolution between the scales A and uj within the EFT. 

For a theory with multiple scales, the procedure is similar. A typical example 
is shown in Fig.|5j which depicts the construction of an EFT at a low scale uj <C 
A2 -C Ai starting from a theory of light fields (j> coupled to heavy fields $1, $2 
(masses of order Ai, A2). One first constructs an EFTi for <f>, $2 (regarded as 
approximately massless) by integrating out the fields $1. This generates a theory 
defined by its coupling constants at the renormalization scale fi ~ Ai. This EFTi 
is then used to RG evolve the couplings down to the threshold /i ~ A2, at which 
point $2 is treated as heavy and removed from the theory. This finally generates 
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RGflow 



\l=A 2 



match 



' RGflow 



(i=CO 



Fig. 5. Construction of a low energy EFT for a theory with scales Ai ^> A2 3> ui. 



an EFT2 for the light fields (f> which can be used to calculate at the scale u>. In 
EFT2, logarithms of u>/ A 2 can be resummed by RG running^. 

Exercise: Calculate the one-loop correction to the Wilson coefficient cg([i). 

1.2.3. Power counting 

In every EFT, there is a rule for determining which operators are needed at a given 
order in the expansion parameter cj/A <C 1. This rule is called power counting. 
It is slightly unusual, in that a typical power counting scheme assigns the same 
counting to tree and loop contributions to a given observable. This is unlike 
the situation in renormalizable field theories where the perturbative expansion is 
equivalent to an expansion in powers of Ti (the loop expansion). 

To see how this works, we develop the power counting rules of our toy EFT. 
In this case ,power counting is simply a matter of keeping track of operator mass 
dimensions^. We simply make the obvious assignment 

P - uj => ~ lo =S> ~ w" 1 . (1.41) 

Since the kinetic operator is leading in the expansion, it should scale as 

J d*x(d,n) 2 ~ (|)° , (1-42) 

2 There are exceptions to the construction described in this paragraph, for example NRQCD, the 
EFT that describes the non-relativistic limit of the strong interactions. Although this theory contains 
several widely separated scales, these are correlated implying that RG running must be performed in 
one stage, rather than in the two-stage picture presented here (23). See 1121 for a pedagogical review. 
This will not be an issue in our discussion of the non-relativistic gravity. 

3 This need not always be the case, as we will see in the next lecture. 
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Fig. 6. Feynman graphs contributing to the -kit — > 7T7t in the effective theory. Graph (a) corresponds 
to one insertion of the dimension eight operator (d^nd^n) 2 . Graphs (b), (c) correspond to insertions 
of a generic dimension ten and twelve operator, respectively. Note that the lj/A power counting 
scheme implies that the one-loop diagram (d) comes in at the same order as the tree graph (c). 



indicating that ir(x) ~ u>, and therefore the scaling of an operator is simply its 
mass dimension. 

Given these power counting rules, we now know how to assign powers of oj/A 
to terms in the expansion of any low energy observable. To see how this works 
let us consider again irir — > itit scattering. By the LSZ reduction formula, the 
S'-matrix element for this process is of the form 



S 



-iki -Xi p\2 



|T7r(ai)--.7r(s 4 )|0). 



(1.43) 



Thus to decide the order in which a given operator contributes to the scatter- 
ing amplitude, we need to know how its contribution to the four-point function 
(T7r(a;i) • • • it{x^)) sales in w/A. 

For example, the leading non-trivial contribution to scattering, from a single 
insertion of the dimension eight operator Og = (d^nd^ir) 2 gives a term whose 
magnitude is 



Fig.Ha) = (T7r(a?i).--7r(a:4) 



g J d*xO s (x) 



A 4 



(1.44) 
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(The subscript means we are calculating the correlation functions on the RHS 
in the free field theory, using Wick's theorem). At higher orders in the expansion, 
the four-point function receives corrections from the diagrams in Fig. H^b), (c), 
(d). In particular, a typical contribution from a dimension ten operator O\o goes 
like 



Fig.Hb) = (T7r(a:i)..-7r(a:4) 



IClQ 

A 6 



d A xO w {x) 



(1.45) 



and likewise a single insertion of a dimension twelve operator On (Fig. |6J C )) 
comes in at relative order w 8 /A 8 . Note however that at this order there is an addi- 
tional contribution from one-loop diagrams with two insertions of the dimension 
eight operator: 



Fig.Hd) = (T7r(a:i)..-7r(a:4) 



A 4 



d 4 x0 8 (x) 



(1.46) 



thus in practice loop graphs can contribute at the same order as tree level inser- 
tions. 



2. Lecture II 

2.1. The binary inspiral as an EFT calculation 

As discussed in the previous lecture, there are three scales involved in the binary 
inspiral problem: the internal structure scale r s , the orbital distance r and the ra- 
diation wavelength r/v. The goal is to calculate gravitational wave observables, 
e.g. the radiated power, arising from physics at the scale r/v. Following our pre- 
vious discussion, a convenient way of doing this is to formulate the problem in 
terms of EFTs. Since the problem has two intermediate scales, the construction 
of the EFT that describes the radiation modes must proceed in stages, as shown 
in Fig. [7] 

The starting point is the theory of an isolated black hole or neutron star cou- 
pled to gravity. Thus in this theory, the relevant degrees of freedom are field 
perturbations (gravitational, electromagnetic) propagating in the BH/NS back- 
ground geometry. Of these degrees of freedom, the field modes of interest in 
the non-relativistic binary dynamics have wavelengths (A ~ r/v) much larger 
than the scale that characterizes the internal structure of the compact object (e.g., 
the Schwarzschild radius r s ). Given our previous discussion, it therefore makes 
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H=l/iv 



pt. particle +GR 



match 



RGflow 



(i=l/r 



bd. state 



mult. + rad. 



match 



RGflow 



(X=v/r - 



Fig. 7. Construction of EFT for binary stars in two stages. 



sense to integrate out the internal structure of the object by matching onto a new 
theory that captures the relevant degrees of freedom. This theory is simply the 
theory of a point particle coupled to the gravitational fields (plus whatever other 
massless fields there are in the problem). We describe how to construct this EFT 
(identifying which modes to keep and constructing the Lagrangian) in sec. 12.21 

This point particle EFT is the correct theory for length scales all the way to 
the orbital radius r. At the orbital scale, the gravitational field can be viewed as 
a superposition of "potential" modes that act over short distances, mediating the 
forces that form the bound system, and "radiation modes" which correspond to 
the gravitational waves that propagate out to the detector. It is convenient once 
again to integrate out the potential modes by matching onto an EFT of radiation 
gravitons coupled to a composite object of size r <C A. This new EFT consists 
of a point particle, together with a set of multipole mass moments generated by 
the mechanical plus gravitational energy of the two-particle bound state. These 
multipoles are coupled to gravity in a way consistent with gauge invariance under 
long wavelength coordinate transformations. We discuss the matching onto, and 
power counting within this theory in sec. 12.41 

Note that in Fig. [7] we have also indicated that the parameters in each of the 
EFTs relevant to the binary exhibit RG flows as a function of a renormalization 
scale. This RG flow, which is purely classical, unfortunately does not arise until 
order v 6 in the binary dynamics, and is beyond the scope of the discussion in 
these lectures. 
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2.2. The EFT for isolated compact objects 

An isolated compact object probed by long wavelength fields can be approxi- 
mated as a point particle. Even if we do not know the internal structure of the 
compact object, it is possible to write down an EFT that describes its interactions 
with external fields which accounts for finite size effects in a systematic fashion. 

Recall from our general discussion that to build an EFT, one needs to do two 
things: 

• Identify the relevant degrees of freedom at the scale of interest. 

• Construct the most general Lagrangian for these degrees of freedom that is 
consistent with the symmetries. 

A black hol^3 interacting with long wavelength gravitational fields can be 
viewed to a first approximation as a point particle probe of the background ge- 
ometry. The degrees of freedom necessary to describe such a system are 

1. The gravitational field g^ix). 

2. The black hole's worldline coordinate x M (A), which is a function of an arbi- 
trary affine parameter A. 

3. An orthonormal frame e^(A) localized on the particle worldline. This de- 
scribes the orientation of the object relative to local inertial frames. It encodes 
how the particle is spinning relative to the gravitational field. 

For the sake of simplicity, I will ignore the effects of spin in the subsequent 
discussion. The construction of effective Lagrangians involving spin degrees of 
freedom can be found in ||251 . 

These degrees of freedom must couple in all ways allowed by the symmetries 
of the problem. These are 

1. General coordinate invariance, x M — > x fl {x). 

2. Worldline reparametrization invariance (RPI), A — > A(A). 

I will make one additional symmetry assumption, again just to keep the discus- 
sion as simple as possible: 

3. SO(3) invariance. This guarantees that the compact object is perfectly spher- 
ical. In particular it has no permanent moments relative to its own rest frame. 

This last assumption, together with the omission of spin degrees of freedom, 
means that the EFT that we are constructing is the appropriate one for describing 
Schwarzschild black holes interacting with external gravitational fields. 



4 Neutron stars may contain additional low frequency modes that must be kept in the point particle 
EFT. See [241 for a review of NS/BH spectroscopy. 
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It is straightforward to write down effective Lagrangians that are invariant 
under these symmetries. To take care of coordinate invariance, we just write 
down Lagrangians that transform as coordinate scalars constructed from g^uix) 
and dx^/dX. A simple way of ensuring RPI is to use the proper time variable 

dr 2 =g liU (x(X))dx' 1 dx v , (2.1) 

as the worldline parameter. Since proper time is physical (i.e., measurable) it 
must be invariant under worldline reparametrizations. 
The effective action consistent with these criteria is then 

s eff[ x *,9nv\ = Seh [g] +Spp[x,g], (2.2) 
where we take the usual action for the gravitational fielci 

S EH = -2m 2 Pl J d 4 x^R(x), (2.3) 

with nip 2 = 2>2itGn, and R(x) is the Ricci scalar. S pp is given by 

S pp = Tti J dr + ■■■ , (2.4) 

where m is the particle mass, and we have suppressed temporarily any possible 
curvature dependent terms in the point particle action. Extremizing S pp gives 
rise to the usual geodesic motion of a test particle in a gravitational field 



=>> 3? + Y^ a RX a x p = x a D a x^ = 0, (2.5) 



(i M = dx^/dr). Including only one power of the curvature, S pp has two addi- 
tional terms 

S pp = — Tft J dr + cji J drR + cy J drR flv x^x 1 ' ■ ■ ■ . (2.6) 

Note however that the leading equations of motion for the gravitational field, 
which follow varying Seh and neglecting sources imply that R flt/ (x) — 0. Thus 
operators constructed from the Ricci curvature are "redundant" operators, and 
can be omitted from S pp without affecting the physical consequences of the the- 
ory. Technically, this is because it is possible to perform a field redefinition of 



5 In addition to the Einstein-Hilbert term, the action for gravity may contain additional powers 
of the curvature suppressed by powers of the scale mpi where new physics is expected to come in. 
These play no role in our discussion. 
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9^v{x) which sets the coefficients of these operators to zero. See the appendix 
for details. 

Since terms with can be omitted, all that is left are operators constructed 
from the Riemann tensor (specifically its traceless part, the Weyl tensor). The 
simplest such operators involve two powers of the curvature and can be written 

as 

S„ = -mJdT + c E J drE^E^ + c B J drB^B^ + ■■■. (2.7) 

The tensors -E M „, B^ v denote the decomposition of the Riemann tensor R^a/3 
into components of electric and magnetic type parity respectively. They are the 
gravitational analog of the decomposition of the electromagnetic field strength 
F^ v into electric and magnetic fields. Explicitly, they are given by 

E„ v = R„ au px a xP, (2.8) 

BfAi/ tfiaPp-E^ R pvX^ . (2.9) 

These tensors are purely spatial in the particle rest frame and for a background 
with = also satisfy E» M = 5% = 0. 

The operators J drE^E^ and J drB^B^ are the first in an infinite series 
of terms that systematically encode the internal structure of the black hole. One 
way to see that these terms describe finite size effects is to calculate their effect 
on the motion of a particle moving in a background field g^. The variation of 
Spp including the terms quadratic in the curvature gives 



dr 



= -S 



c E J drE^E^ +c B j drB^B^ + ■ 
x a D a x tl ^ 0. (2.10) 



In other words, due to the curvature couplings, the particle no longer moves on a 
geodesic. However, geodesic deviation implies stretching by tidal forces, which 
occurs when one considers the motion of extended objects in a gravitational field. 
In fact, an explicit matching calculation, described in Sec. 12.61 predicts that the 
coefficients ce,b ~ m 'pi r s^ vanishing rapidly as the size of the black hole goes 
to zero. 

The EFT in Eq. ( 12.21 ) describes the dynamics of one extended object in a grav- 
itational field. To describe the motion of several objects, we simply include a 
separate point particle action for each 

Seff[x a ,g] = Seh [.<?] + ^2Sp p [x a ,g], (2.11) 

a 

where the index a runs over all the particles moving in the field g^ u . 
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2.3. Calculating observables 

In principle, Eq. ( 12.1 It correctly captures the physics of an arbitrary system of 
extended objects for scales fi < l/r s . It can therefore be used to calculate all 
observables measured by gravitational wave detectors. 
If we decompose the typical gravitational wave signal as 

h(t)= A(t) cos <f>(t), (2.12) 

then interferometric detectors, such as LIGO/VIRGO and LISA, are particularly 
sensitive to the phase <fr(t) ~ 2 J* oItuj{t) of the gravitational wave. This is usu- 
ally calculated using the "adiabatic approximation": Consider a non-relativistic 
binary inspiral, and suppose we have calculated to some order in the v <C 1 
expansion the quantities 

E(v) = mechanical energy of binary. 

P{v) = power emitted in gravitational waves, 

as functions of the orbital parameters, such as the velocity v(t), or equivalently 
the orbital frequency uj(t). Energy conservation 

then gives a differential equation that can be used to solve for u(t) which in turn 
gives the frequency of the GW wave signal uiQyy(t) ~ 2<j(t) and consequently 
the phase <p(t) (the factor of two arises from the helicity-two nature of the gravi- 
ton). We implicitly did this calculation to leading order in v in Sec. 11.11 There 
we found 



1 1 



v{tf v{Uf 



(2.14) 



for equal mass stars in a circular orbit. 

Although the observables we are interested in calculating are purely classical, 
let us pretend for the moment that we are doing quantum field theory. Write 

mpi 

and calculate the functional S e ff(x a ) defined by the path integral 

exp [iS e f f (x a )] = J Dh fJlU (x) exp [iS EH (h) + iS pp (h, x a )\ , (2.16) 
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with the particle worldlines x%(t) that source /i M „ held fixed. It turns out that 
S e ff(x a ) is a generating function for the quantities of interest. In particular, the 
classical limit of the variation 



6\ReS eff (x a )]=0 



(2.17) 



gives rise to the equations of motion for the worldlines x£(t). From this one can 
derive an energy function E(v) in the usual way. Furthermore 



-lmS eff (x a ) = - 



dEdVL- 



d 2 T 
'dfldE' 



(2.18) 



measures the total number of gravitons emitted as worldlines (r) evolve over 
a time T — > oo. Although graviton number is not a well defined observable 
classically, the power emitted can be obtained by integrating the differential rate 
dT over the energy of the emitted graviton 



P = tot. power 



dEdn e- 



d 2 T 
dVLdE ' 



(2.19) 



To see how the calculation of observables via Eq. ( 12.161 ) works, we will con- 
sider a toy gravity model consisting of a scalar graviton field 4>(x) interacting 
with several point particles. The scalar <f> couples to the point particles with 
strength proportional to mass 



S = \ ( d^xd^d^cj) -^2m a f dr a 

J a J 

or equivalently 

dr„. ■ 



1 



S = - ^ m a 

a 

with 

j(x)^-J2 



d*x 



d^(j)+J{x)(j){x) 



dT a S 4 (x - x a ). 



(2.20) 



(2.21) 



(2.22) 



2\f2mpi 

Then the functional that generates the observables in this model is given by 

S e ff(x a ) = ~y^m a J dr a - i\nZ[J], (2.23) 



where 



Z[J} = J D(j)(x)exp i J dtxha^d^ct) + J(x)<f>(x) 



(2.24) 
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This is a Gaussian integral so it can be easily calculated explicitly. Up to an 
irrelevant constant, the result is 

In Z[J] = -\f d i xd i yJ{x)D F {x - y)J{y). (2.25) 

In the case of real gravity, with non-linear self interactions, the analog of Z[J] 
will not have such a simple expression. Nevertheless, the perturbative expansion 
of the generating function in real gravity has a simple diagrammatic interpreta- 
tion. As a warmup to constructing the diagrammatic rules in real gravity, it is 
useful to recall how Eq. (12.25b is recovered by summing up Feynman diagrams. 
Treat the coupling J d 4 xJ(x)(f>(x) perturbatively, and introduce the Feynman 
rule 

~T~ - i /^-"-=- i ?i^/^-" 1 -- 

(2.26) 

Then by expanding the interaction term in Eq. d2.241 >. we see that Z[J] has a 
diagrammatic expansion in terms of "ladder diagrams" 



™ - ~T^TT+TTT 



I 



exp S . (2.27) 



Note that the intermediate particle lines have no propagators associated with them 
(they just depict the time evolution of the particle worldlines), so diagrams with 
multiple "rungs" are simply products of the diagram with a single scalar ex- 
change. Specifically, the diagram with n intermediate graviton propagators is 
given by 



JZZT-h. 



d 4 xd 4 y(iJ(x))D F (x-y)(iJ(y)) 



(2.28) 



where the factor of l/2 n n! appearing here is the symmetry factor associated with 
the diagram. Summing up contributions from diagrams with any number n of 
scalar lines, we see that the series in Eq. ( 12.271 ) reproduces Eq. ( 12.251 ). 

Notice that the exponentiation of the diagrams contributing to Z[J] implies 
that S e ff(x a ) receives a contribution only from diagrams that remain connected 
after the particle worldlines are stripped off, which in this theory is just the di- 
agram with a single graviton exchange. The fact that 5 e //(x a ) is given by di- 
agrams that do not break up after removing all worldlines remains true in real 
gravity, where the graphs can have graviton self-interaction vertices. 
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The effective action in the toy gravity model is thus given by 

S e ff(x a ) = ~y^m a J d,T a + ^ T^T* y dT a dnD F (x a - x b ). 



O . & 



(2.29) 



Eq. ( 12.291 ) contains all the information about the classical particle dynamics. For 
illustration, let's evaluate it for the case where the motion is non-relativistic. In 
this case the integral 



dT a dT b D F (x a - Xb) = / dradn 



d A k 



(2tt) 4 k 2 + it 



-ik- {x a —Xb) 



(2.30) 



can be split up into contributions from two regions of momenta: 

• Potential: This is the region corresponding to gravitons with spacelike mo- 
menta of the form 

fc/* _> (fc° „ _,k~ -), (2.31) 

r r 

where r is the typical separation between the particles and v <C 1 is the typical 
three-velocity. Gravitons with this momentum configuration mediate nearly in- 
stantaneous exchanges between the point particles. Note that potential gravitons 
can never be on-shell, k 2 ^ 0. Therefore they never contribute to the imaginary 
part of Eq. ( 12.29b . which happens when the it term in the propagator becomes 
important, i.e., when k 2 = 0. 

• Radiation: Radiation gravitons correspond to on-shell modes with 



k" -> (fc° 



,k~-). (2.32) 

r 

Because these modes satisfy k 2 ~ 0, they are responsible for giving rise to the 
imaginary part of S e ff (x a ) and therefore the radiation that propagates out to the 
detector. 

Since only potential modes contribute to the real part of S e ff(x a ), we may 
calculate in the NR limit by expanding the propagator as 



1 



k 2 



Using 



1 

k 2 " 



k 2 

1+ *0 
k 2 



^[1 + ^ 2 )] 



(2.33) 



d 4 k 

(27T)4 



-ik-x 



l 

47r|x 



■S(x°), 



(2.34) 
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and expanding dr a — dx°^/l — v 2 ~ dx° [l — , we find to leading order 



in v 



*eS eff K] = \ E / dx ° m °< - \ E / ^ ^ff + • • • , (2.35) 

L a J Z a,b J 1 a fc| 

which is just the Lagrangian for classical NR particles interacting through a New- 
tonian potential. Note that the second term contains divergent self-energy con- 
tributions whenever a = b in the sum. These divergences can be absorbed by 
renormalization into the particle masses. Formally, they can just be set to zero, 
by evaluating the momentum integral in dimensional regularization using the for- 
mula 



d d k 



— tk-: 



i 



1 T(d/2 - a) 



(k 2 )« (4ir) d / 2 r(a) 



2\ a ~ d / 2 



(2.36) 



and taking the limit x — > before setting d = 3. 

To calculate the imaginary part of S e ff(x a ) we use 



Im 



fc 2 



-inS(k 2 ), 



(2.37) 



which physically has the effect of ensuring that only on-shell particles contribute 
to the radiated power. Then 



lmS ef f(x a ) 



1 



/ (S^i E m «/ dr « e "* Xo 



16m£, y (2tt) 3 2|k| 



(2.38) 



fc°=|k| 



and the differential power is 
dP 



= 1 Qiir 2 

dttd\k\ ~ Tin 2 



(2.39) 



fc°=|k| 



This is exactly what one would find classically, by solving the wave equation for 
4> with source term given by Eq. (12.221 1 and then calculating its energy-momentum 
tensor at an asymptotically large distance from the source. It can also be calcu- 
lated from the tree-level amplitude for single graviton emission from the source 
particles 



*E 



2mpiy/2 



dr a e l 



(2.40) 
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The differential probability for the emission of one graviton is then 

d 3 k 1 



<iProb(vac. 



= \A(yac. <ji)Y 



(2tt) 3 2|k| 



(2.41) 



Dividing this result by the observation time T — ► oo to get a rate, and multi- 
plying by the graviton energy |k| to convert to power reproduces Eq. (12.391 . So 
indeed the functional integral formula for S e ff{x a ) knows about the classical 
observables for the particle ensemble. 

Exercise: In the case of electrodynamics coupled to point particles, 



S =-\J d^xF^F^ -J2 m « J d T a + J2 eQa J d < A ^ (2 



42) 



show by calculating lmS e f / (x a ) that the radiated power in photons is given by 

2 



dP 



1 a 



dCld\k\ T4ir 2 



(2.43) 



feo=|k| 



where a — e 2 / (4tt). Verify that, upon time averaging, this reproduces the usual 
electric dipole radiation formula in the non-relativistic limit. 

Gravitational radiation in general relativity can be calculated by a straight- 
forward generalization of the methods used for the scalar gravity model. Real 
gravity is of course more complex, but the complications are mainly computa- 
tional. If we expand as in Eq. ( 12.15b and plug into the gravitational action 
we find an infinite series of terms that are schematically of the form 



2mp ; J d i x y /gR(a 



d 4 x 



{dhf 



h(dh) 2 h 2 (dh) 2 



mpi 




(2.44) 



leading to graviton self-interactions with Feynman vertices containing any num- 
ber of graviton lines. We will not need the detailed form of these terms in what 
follows. A derivation of the Feynman rules for gravity can be found in the lec- 
tures by Veltman ||26l . See also [19|. In addition, the gravitational field has 
non-linear interactions with the point particle, e.g., 
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(a) (b) (c) 

Fig. 8. The first few diagrams contributing to S e ff(x a ) in Lorentz covariant perturbation theory. 



8m 2 pl 



df ( h^xPx") 2 + ■ ■ ■ , (2.45) 



where dr 2 = r\ jiv dx ii dx v . This gives rise to vertices 



■\f\j 



s 3 

h 



(2.46) 



As a result of this, the diagrammatic expansion of S e ff(x a ) in the gravita- 
tional case is much richer than in the simple model considered above. The first 
few terms in the expansion are shown in Fig. [S] The problem with this (covari- 
ant) form of the perturbative series is that it is not optimal for the calculation of 
observables in the v -C 1 limit. For example, consider the diagram in Fig.[8fc). 
At what order in v does it contribute? In the language of the previous lecture, 
we do not have a velocity power counting scheme for the diagrams appearing in 
Fig. [8] Since in the end we are interested in computing gravitational wave signals 
to a fixed order in v, it is important to develop a set of rules that assigns a unique 
power of v to each diagram in the theory. 



2.4. Integrating out the orbital scale 

The reason why the diagrams do not scale as definite powers of v is that there are 
multiple scales in the theory. Although we have successfully integrated out the 
internal structure scale r s , the momentum integrals in the diagrams of Fig.[8]have 
contributions both from potential (v/r, 1/r) and radiation (v/r, v/r) gravitons. 
But potential modes are never on-shell (they have k 2 ~ 1/r) and do not belong 
in the effective theory at long distance scales. 

To remedy this problem, we will explicitly treat the potential and radiation 
modes separately, by splitting up the gravitational field as 



(2.47) 
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Here, the field represents the long-wavelength radiation modes. Schemati- 
cally, it satisfies 

v - 

dah^v ~ -V" ( 2 - 48 ) 
r 

meaning that the field varies slowly over spacetime, with a typical length 
scale r/v. It can therefore be regarded as a slowly varying background field in 
which the potential modes H^ v , with 

doH^ ~ diH^ ~ -H^ u , (2.49) 

propagate. Actually it is useful to re-write in terms of its Fourier transform 

/<i 3 k 
—-^e^H^ixo). (2.50) 

This allows us to explicitly disentangle hard momenta 

k~l/r, (2.51) 
from the long wavelength scale 



T 



(2.52) 



and to treat all derivatives acting on fields on the same footing, ~ v/r. 

To derive an EFT that has manifest velocity power counting rules, we need to 
integrate out the potential modes. The basic idea is the following: we calculate 
S e ff(x a ) (and hence all observables), in two stages. In the first stage, we perform 
the integral over Hy^ v , 



e iS N R.(h,x a ) _ 



= J Z?ff k ^(xV SBH(S+H)+Spp( * +H ' ;i;a) (2-53) 

where the field is treated as a background field. This gives the quantity 
SnrQi, x a ), which formally contains the two-body forces between the point par- 
ticles, written as an explicit expansion in powers of v, and the couplings of the 
particles to radiation. In this theory, the short distance scale r appears explicitly, 
in the coefficients of operators. In other words, Eq. ( I2.531 l is simply a formal way 
of doing the matching to the long-distance EFT containing h^ v and the particle 
worldlines. Given SNn(h, x), the second stage in the calculation is to perform 
the functional integral over h^. Feynman diagrams are easier to calculate in this 
theory than in the full theory, as momentum integrals only contain the scale r/v. 

Note that in Eq. ( I2.531 l we have dropped gauge fixing terms necessary to make 
sense of the path integral. We will not need the explicit form of such terms in the 
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subsequent discussion, but it is perhaps worth noting that it is very convenient to 
do the gauge fixing in a way that preserves the invariance under diffeomorphisms 
that transform the background metric h^ u (this is called the background field 
method 11271 ). If such a gauge fixing scheme is chosen, the action for h^ v is 
guaranteed to be gauge invariant. This places strong restrictions on the form of 
the EFT that describes radiation. 

Diagrammatically, SmrQi, x a ) is given by a sum over diagrams that have the 
following topological properties: 

1 . Diagrams must remain connected if the particle worldlines are stripped off. 

2. Diagrams may only contain internal lines corresponding to propagators for 
the potential modes Hi if2[/ (xo). Diagrams cannot contain external potential 
graviton lines. 

3. Diagrams may only contain external lines. Since the functional integral 
in Eq. ( 12.531 ) views the field h^ v as a background, diagrams cannot contain 
propagators corresponding to internal radiation graviton lines. 

The point of splitting the original graviton into the new modes h^, H-^^ 
is that the diagrams in terms of these new variables can be assigned definite pow- 
ers of the expansion parameter v. The power counting rules for determining how 
many powers of v to assign to a given diagram follow simply from the fact that the 
three momentum of a potential graviton scales as k ~ 1 /r, since this is the range 
of the force it mediates, and that the spacetime variation of a radiation graviton 
is x^ ~ r/v. These two observations are sufficient to consistently assign powers 
of v to any term in the action, and by extension to the Feynman rules. 

First, let's determine the power counting for diagrams that only contain po- 
tential modes. We can fix the scaling of H^ v in terms of kinematic variables by 
looking at its propagator. To obtain the propagator, we need the quadratic terms 
in H^ u from the gravitational Lagrangian. Including an appropriate gauge fix- 
ing term, whose form is not necessary for our purposes here (see [7 1 for details) 
the relevant terms are 



1 f , n d 3 k 



S H * =-\f dx° 



(2tt) e 



k 2 

k Hkp U H_k — —HkH-k 



(2.54) 



where = H a ak. Terms with time derivatives have been dropped, as do is 
down relative to k by a power of v. Because of this, the potential propagator is 
instantaneous, 

(THk^ix^H^O)) = -±(2n) 3 S 3 (k + q)5(x a )P^, af3 , (2.55) 
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where the tensor structure is given by 



(2.56) 



Eq. (12.55b is essentially the Fourier transform of the Newton potential, as one 
would expect. Since P^-ap is 0(1), the scaling of is simply 



Hi 



so that 



4 

= r v, 



(2.57) 



(2.58) 



Given the scaling rule for H^ v , it is now possible to assign powers of velocity 
to diagrams with no external radiation gravitons. For example, the coupling of 
i/koo t° a NR particle is 



2mp( 



j ^ k ik-x(i°) tj 



m ( r 
XI — ] 

mpi \v. 

m 



x I rV' 2 



V 



-1/2. 



mpi 



The scaling of m/mpi is fixed by the virial theorem, 



mv r = Lv, 



(2.59) 



(2.60) 



' pi 



where we have used Gat ~ l/m Pl , and have introduced the orbital angular 
momentum L = mvr. The interaction in Eq. ( 12.591 ) then scales as 



2mpi 



(2tt)3 



kOO 



(2.61) 



so the exchange diagram with two insertions of this term, which leads to the 
Newton potential between the point particles, scales as a power of the angular 
momentum L, 



(2.62) 
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It can be shown that any diagram without external radiation lines scales like 
L n v m , with n < 1 and in > 0. The bound on n is saturated by diagrams 
without graviton loops. Thus L is the loop counting parameter of this theory, and 
since we are interested in the classical limit ?l — > 0, it is a good approximation to 
drop diagrams with graviton loops (formally, such diagrams are down by powers 
of h/L < 1). 

Power counting rules for potential graviton self-interactions can be derived 
along similar lines. Consider the potential three-graviton vertex. Expanding the 
gravitational Lagrangian to cubic order in H^^, one finds terms with the struc- 
ture 



S H 3 



1 



mpi 



dx°(2Tr) 3 5 3 {J2 k k; 



3 ^3 



(2.63) 



where we have inserted the expansion Eq. ( 12.501 ) into the Einstein-Hilbert action 
and performed the d 3 x integral . The factor of k 2 comes from the fact that every 
term in the action has two derivatives on the metric. From this we find 



J- x n x / 1 

mpi \v 



1 



and therefore the diagram 



(2.64) 



VIY x x VI = Lv 2 , 



(2.65) 



gives rise to a term in the two-body potential that is suppressed by v 2 relative to 
the leading order Newtonian exchange diagram. The full set of potentials at order 
v 2 arise from the diagrams in Fig. [9] In addition to the three graviton term just 
discussed, there are additional velocity suppressed vertices from the expansion 
of the particle proper time Lagrangian in powers of v, 



S i, 



mpi 



dx° 



-■^h 00 - h 0i Vi - ift-oov 



2 hjji 

2 



(2.66) 



The sum of these diagrams gives the corrections to the two-body NR Lagrangian 
at order v 2 



1 
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|xi - x 2 | 
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Fig. 9. Diagrams contributing to the order v 2 corrections to the two-body potential. 



(vi • xi 2 )(v 2 ■ x i2 ) 



Xi - x 2 



G%m 1 m 2 (m 1 + m 2 ) 
2|x 1 -x 2 p • 



(2.67) 



where X12 = Xj — x 2 . This Lagrangian was first obtained by Einstein, Infeld and 
Hoffman l28l in 1938. In this equation, the first term is just the first relativistic 
correction to the particle kinetic energies, the second term arises from diagrams 
with a single graviton exchanged by velocity-dependent vertices, and the last 
term is from the diagrams of Fig. [3c), (d). Diagrams at higher orders in v can be 
calculated and power counted by the same methods outlined here. 



2.5. Radiation 

Integrating out the non-dynamical potential modes generates the gravitational 
forces between the non-relativistic particles. These arise from diagrams with no 
external factors of the radiation field h^. The functional integral in Eq. (12.531 1 
also generates, from the diagrams with one or more external radiation graviton, 
the couplings of matter to radiation. 

Incorporating radiation is fairly straightforward. First note that the propagator 
for the radiation field, in a suitable gauge, is given by 

(TV(x)^(O)) = j ^yk^+ie e ~ tk ' XP ^- (2M) 
Since fc M ~ v/r, Eq. ( |2.68t implies that radiation modes should scale as 

V ~ -. (2.69) 

r 

This rule allows one to power count terms in the action containing the radiation 
field. However, in order to obtain an EFT for radiation that has manifest velocity 
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power counting, the decomposition of the graviton into potential and radiation 
modes is not sufficient. It is necessary also to multipole expand the couplings of 
the radiation field to either potentials or to the point particles, 

V(x,z°) = V(X,s°) + <5x < a i V(X,s ) + i*xW'/i /il/ (X,a; ) 

+ • • • , " (2.70) 

where X is an arbitrary point, for example the center of mass of the multi -particle 
system, X cm = ^ a m a x a / m a- The reason for this is familiar in the case 
of couplings to matter sources: consider the amplitude for graviton emission by 
an ensemble of point particles (v% is the four- velocity, e flu (k) is the graviton 
polarization) 

iA = T , r£ =~iT-£±- f dr a e ik ^v^;M- (2-71) 

a a 

The final state graviton is on-shell, so k° = | k | . On the other hand, if the particles 
are non-relativistic, x° ^> |x a | (measuring positions relative to the center of 
mass). Then 



l + ik-x a -i(k-x Q ) 2 + . 



(2.72) 



Note that if |k| ~ v/r and |x a | ~ r, k • x Q ~ v. So the exponential contains 
an infinite number of powers of v, and the amplitude does not scale homoge- 
neously in velocity. This is avoided by multipole expanding at the level of 
the Lagrangian. Likewise, if we consider a graph containing potential as well as 
radiation modes 

Fig.[Ha~ - * = -l[l-2p ■!£+■■■], (2.73) 

(p + ky p z 

we find that because |p| ~ 1/r and |k| ~ v/r, the propagator for the potential 
mode with three-momentum p + k contains an infinite number of powers of 
v. To ensure that this does not happen, it is necessary to arrange that radiation 
gravitons do not impart momentum to the potential modes. This is exactly what 
occurs if we plug in the expansion Eq. ( 12.701 ) into the couplings generated by the 
gravitational Lagrangian. The necessity to multipole expand potential-radiation 
couplings was first pointed out in the context of non-relativistic gauge theories in 
ref. 1291. 

Given these ingredients, we are finally in a position to calculate the matching 
coefficients for the radiation graviton EFT at the scale fi ~ 1 jr. Matching boils 
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(p°,p) (p° + k°, P + k) 

Fig. 10. The interaction of a potential and a radiation mode. 

down to calculating Eq. (I2.531 l. or equivalently comparing diagrams with any 
number of external radiation modes in the full theory and the EFT. For example, 
consider the graphs with one external h^. In the theory with both potentials and 
radiation, the leading one is 



E 



±-fdx°J2™ahoo(0,x ). (2.74) 



2m p i 



In this equation, the double lines on the graph on the right indicates that this is 
a vertex in the EFT below /i ~ 1/r, where the separation r between the point 
particles (a short distance scale) cannot be resolved. It is the vertex obtained by 
matching at the potential scale, and it corresponds to the term in the EFT shown 
on the right hand side of the equation. At this order, the matching is just the ze- 
roth order multipole expansion (performed about the center of mass), i.e. which 
generates the coupling of hoo to the mass monopole m = J2 a m a- At the next 
order, the effective theory vertex contains terms both from the multipole expan- 
sion and from the explicit velocity dependence of the point particle gravitational 
couplings. The result is 



E 



2m p i 



dx° [X 



,, m dih, 



oo 



2P% t ] , (2.75) 



where P = m a'^-a is the total linear momentum of the system. In the CM 
frame X cm = P cm = 0, so this coupling vanishes, i.e., there is no dipole radia- 
tion in general relativity. 

Things are more interesting in the one-graviton sector at order v 2 . The graphs 
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at this order are 



(2.76) 



where the first graph depicts the terms from the multipole expansion and from 
the point-particle couplings containing two powers of the velocity, the second 
graph (an almost identical mirror image graph not shown) comes from the two- 
graviton couplings of the point particle, and the vertex in the third term is from 
the hH 2 terms in the Einstein-Hilbert Lagrangian. The sum of these graphs gives 
the second order terms in the EFT, S v z [h] — J dx°L v 2 [h], where 



L V 2 [h] 



2mpi 



— ^2 ?71aX^X J Q i? 0i C 



|xi - x 2 | 



1 



2m 



Pi 



-eijklikdjhoi 



(2.77) 



Here, the first term is just the coupling of hoo to the Newtonian energy of the 
two-particle system. This term can be regarded as a kinetic plus gravitational 
correction to the mass monopole of the source 



1 r 



|xi - x 2 | ' 



(2.78) 



The second term is a coupling of the graviton to the total mechanical angular 
momentum of the system, L = ^ a x a x m a v Q . Both the mass monopole and 
the angular momentum are conserved at this order in the velocity expansion, 
and therefore cannot source on-shell radiation. Rather they source the long range 
static gravitational field of the two-body system. The last term, on the other hand, 
is a coupling of the moment m a x a x a to tne (linearized) Riemann tensor of 



the radiation field. Note that for on-shell radiation i?oo 
traceless moment 



Rn 



E 



0, so only the 



(2.79) 



is a source of gravitational waves. It is interesting to note that in order to ob- 
tain Eq. ( I2.771 i. we needed to include not only the multipole expansion of the 
couplings of h^ u to the point particles, but also diagrams with graviton self- 
interactions. This is forced upon us by the power counting rules of the theory. 
It would be simply inconsistent to drop the second two graphs in Eq. ( 12.771 and 
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keep only the first. Thus even the leading order quadrupole coupling to radiation 
is sensitive to the non-linearities present in relativistic gravity. Physically, this 
is perhaps not too surprising, as the graviton couples to all sources of energy- 
momentum, including the energy-momentum stored in the gravitational field of 
the particles themselves. 

Exercise: Use the power counting rules to check that 




(2.80) 



so that the second two graphs in Eq. ( I2.761 l each scale as \fLv 5 / 2 . This is sup- 
pressed by v 2 relative to the graph in Eq. ( 12.741 ). 

Working out higher order terms in the EFT is similar, and there is no concep- 
tual problem in carrying out the expansion to any desired order in v. Once we 
have calculated the terms in Snr[Ii], it is straightforward to compute physical 
observables such as the gravitational radiation power output. As discussed in 
Sec. 12. 31 these can be obtained from 

exp[z5 e// (x a )] = J Dh^(x)e lS -^\ (2.81) 

which is simply the sum over diagrams, calculated using the Feynman rules from 
SNFi[h], that have no external graviton lines. For example, the leading order 
quadrupole radiation formula follows from 

Im C^^Z, => P = Q./lj ■ (2.82) 

where the brackets denote a time average. In order to compute the power, one 
needs to know the time evolution of the moment Qij(t). This can be obtained 
by solving the equations of motion for the x a (i), which follow from the terms in 
Snr[7i] with no powers of h^. 

2.6. Finite size effects 

At some order in the velocity expansion, the internal structure of the binary star 
constituents play a role in the dynamics. In the case of black hole binaries, it is 
easy to determine the order in v for which this happens. As discussed in Sec. 12.21 
finite size effects are encoded in non-minimal couplings to the curvature. Such 
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terms are built out of invariants constructed from the Riemann tensor, the simplest 
ones being those of Eq. d2.71 >. 

S = ■ ■ • + c E [ drE^E^ + c B [ drB^B^ + ■■■. (2.83) 



The coefficients ce.b can be obtained by a matching calculation. One simply 
compares some observable in the point particle theory to the analogous quantity 
in the full theory, adjusting the coefficients to ensure that the two calculations 
agree. 

For the purposes of calculating the coefficients, ce,b, a convenient observable 
is the S'-matrix element for a graviton to scatter off the BH background. In the 
EFT, this amplitude is schematically of the form 



i^-u; 4 + ---. (2.84) 

m pi 



Here, to is the energy of the incoming graviton. For the point particle EFT to be 
valid, it must satisfy r s u <C 1. The factor of 1 /m Pl is due to the fact that we are 
looking at a two graviton process. We have not shown, for example, the leading 
order term in the amplitude, due to graviton scattering off the field created by the 
mass term — m J dr in the point particle action. Although the leading effect of 
ce,b is from interference with this term, the above equation predicts that the total 
scattering cross section contains a term 



c 2 

o-(ui)eft + -^^ 8 + ---. (2.85) 



The only scale in the full theory is the BH radius r s , so the cross section must be 
of the form 

a{w) BH =r 2 J(r s u;), (2.86) 

where f(r s u>) is a function that for r s oj 1 can be expanded in powers of r s oj 
(possibly times logs). Thus we expect to find a term 

o-(u;)bh + rl°iu s +■■■, (2.87) 

and therefore ce,b ~ m p; r s ■ Using this scaling, the power counting rules in the 
non-relativistic limit indicate that the finite size operators first contribute to the 
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dynamics through their effect on the two-body forces, through the diagram 



so is down by ten powers of v relative to the Newton potential. One concludes 
that finite size effects are completely irrelevant for binary inspirals. Actually, this 
is not completely true, as realistic black holes usually have non-zero spin, and the 
inclusion of spin tends to enhance finite size effects [25 1 . Also, absorption by the 
black hole horizon, which necessitates the introduction of additional worldline 
modes in the EFT [ 30 1 arises at order v s for non-spinning black holes [31] and at 
order v 5 [32] when the black hole has spin. 

3. Conclusions 

These lectures presented an introduction to the basic ideas of EFTs and their use 
in understanding the evolution of coalescing binary stars in the non-relativistic 
regime. In treating such systems, one encounters a wide range of physically 
important length scales, from the internal size of the binary constituents to the 
gravitational radiation wavelength. In order to make sense of physics at all these 
scales, it is extremely convenient to formulate the problem one scale at a time, 
by constructing a tower of EFTs as outlined in these lectures. 

One topic that was not discussed in these lectures is the issue of ultraviolet 
divergences in the non-relativistic expansion. Such divergences arise even clas- 
sically, in the computation of the Feynman diagram that contribute to a given 
observable. Physically, the presence of these divergences can be attributed to the 
singular nature of the point particle limit. In the EFT approach, these can be 
handled by the usual regularization and renormalization procedure found in text- 
books: divergences get cutoff using any convenient regulator (e.g. dimensional 
regularization) and renormalized into the coefficients of operators in the point 
particle Lagrangian. Since by construction the point particle EFT has all possi- 
ble operators consistent with symmetries, all short distance singularities can be 
removed from the theory. At high enough orders in the expansion, some of the 
EFT coefficients are logarithmically renormalized, leading to RG flows that can 
be exploited to compute terms in the velocity expansion that are logarithmically 
enhanced. See [7| for more details. 

We expect that the technology of EFTs should be applicable to other astro- 
physical sources of relevance to gravitational wave physics. For example, it is 
expected that LISA will detect gravitational waves from the motion of a small 



CE,B 




(2.88) 
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object, i.e. a neutron star, around a super-massive black hole. In this case, there 
is a small expansion parameter, namely the ratio of the neutron star size to the 
curvature length of the black hole background. It seems natural to construct the 
expansion in this parameter using EFT methods. 

The point-particle EFTs discussed here may also potentially play a role in the 
problem of tracking the evolution of colliding black holes in numerical relativ- 
ity. In particular, the EFT approach to parameterizing black hole internal struc- 
ture could be useful in addressing the question of how to systematically handle 
the curvature singularities present in numerical simulations. More generally, the 
methods presented in these notes are useful for understanding the dynamics of 
extended objects interacting with long wavelength fields. This sort of situation 
arises in both formal and phenomenological applications. Work on the ideas out- 
line here is underway. 
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Appendix A. Redundant operators 

In general, an effective Lagrangian contains all operators constructed from the 
light degrees of freedom that are invariant under the symmetries of the low en- 
ergy theory. In practice, some operators are redundant, and may be dropped 
without altering the physical predictions of the EFT. In fact, operators that vanish 
by the leading order equations of motion may be omitted from the list of opera- 
tors appearing in the effective Lagrangian. Such operators are sometimes called 
redundant operators. 

The basic reason for this is simple. If an operator vanishes by the leading 
order equations of motion (i.e., it vanishes "on-shell"), this means that one can 
redefine the EFT fields in such a way that the Lagrangian written in terms of the 
new fields does not contain the redundant operator. But field redefinitions have 
no effect on physical observables (S'-matrix elements) so the Lagrangian without 
the redundant operator is equivalent to the original Lagrangian. 
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To see how this works in practice, consider in the EFT for the massless tt field 
discussed in Lecture I, the dimension-ten operator 

O w (x) = (d\)F[w,d li w}, (A. 1) 

with 

F = {d a dpTT)d a Tid P TT. (A. 2) 

0io potentially contributes to to the amplitude for nn — > nn scattering at order 
(w/A) 6 . Because to leading order in the uj/A expansion the equation of motion 
for tt is d 2 Tr(x) = 0, this operator actually gives a vanishing contribution to 
scattering. This is obvious at tree level, but in fact Oio(x) may also be dropped 
from loop diagram contributions to scattering amplitudes as well. 
Starting from the Lagrangian including Oio(x) 

S=\f d 4 x(d^) 2 + --- + ^ J d 4 xO(x) + ---, (A3) 

define a new field tt (x) by 

n(x) = jf(x) + Sn(x). (A. 4) 

Then 

J d A x(d^f = j d 4 x(d„Tf) 2 -2 J d 4 xd 2 7fSTr(x) + 0(Stt 2 ). (A. 5) 
So if we choose 

Sir(x) = ^F[W,d ll W], (A. 6) 

we see that the shift in the leading order (dimension four) term in the Lagrangian 
precisely cancels the term with the 0\q(x) in the original Lagrangian. Note how- 
ever that the effects of this operator are not completely gone. For example, insert- 
ing the field redefinition into the dimension eight operator Os(x) = {d^d^ir) 2 
generates operators of dimension fourteen and higher 

(d^d^f -> (d^d^f + ^ (d„7Td»?r)d v ?rd»F +■■■ (A. 7) 

However, when constructing the original EFT, one must include all possible oper- 
ators consistent with the symmetries at every order in the w/ A expansion. Thus, 
the shift in Os(x) due to the field redefinition can be absorbed into the coeffi- 
cients of a operators that are already present in the theory. The same is true for 
the substitution of the field redefinition into any other operator, and we conclude 
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that in terms of the new variables, the effects of the redundant operator 0\q(x) 
are completely spurious. 

The fact that operators that vanish on-shell may be omitted is very general and 
very useful. A systematic discussion is given by Politzer in [ 33 1 . See also ll34l . 
We used this result in the EFT that describes an extended object coupled to grav- 
ity. The effective theory for the point particle contains the terms 

cr / drR + cv / drR^x^x", (A. 8) 



which are redundant due to the fact that at leading order (ignoring coupling to 
sources) the Einstein equations imply R lllJ (x) = 0. Suppose we re-define the 
metric appearing in the original Einstein-Hilbert Lagrangian, — > g M „ + 5g^ v , 
with 



Sg^(x) = - 2 / dr 



1 f 6 4 (x-x(t)) 



2m pi J V9 



'{£,R~\iv)g^+S,v±^x v 



.(A. 9) 



When plugged into the gravitational action Seh this redefinition induces the 
shifts 

cr.v -> c R y + £, R , V . (A. 10) 

Thus by adjusting we can make the coefficients cr y whatever we like, 
including zero, without affecting the physical predictions of the theory. 
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